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Simulations of scanning tunneling microscopy measurements for molecules on surfaces are tradi- 
tionally based on a perturbative approach, most typically employing the Tersoff-Hamann method. 
This assumes that the STM tip is far from the sample so that the two do not interact with each other. 
However, when the tip gets close to the molecule to perform measurements, the electrostatic inter- 
play between the tip and substrate may generate non-trivial potential distribution, charge transfer 
and forces, all of which may alter the electronic and physical structure of the molecule. These effects 
are investigated with the ah initio quantum transport code smeagol, combining non-equilibrium 
Green's functions formalism with density functional theory. In particular, we investigate alkanethiol 
molecules terminated with either CH3 or CF3 end-groups on gold surfaces, for which recent experi- 
mental data are available. We discuss the effects connected to the interaction between the STM tip 
and the molecule, as well as the asymmetric charge transfer between the molecule and the electrodes. 



I. INTRODUCTION 

Scanning tunnelling microscopy (STM)^ is an invalu- 
able surface characterization tool with multiple applica- 
tions in molecular electronics. Typically, an atomically 
sharp probing tip constructed from heavy metals such 
as platinum, tungsten or iridium scans a substrate by 
means of a tiny tunneling current. The method can 
be used as a topographic tool to map the positions of 
atoms, defects and molecules on surfaces, as a spectro- 
scopic tool to probe their local density of states, and 
more recently in its spin-polarized version as an ultra 
sensitive magnetometer—. Furthermore STM can also be 
used as an atomic fabrication tool for depositing atoms 
and molecules on a surface so as to form nanoscale 
devices^. Additionally, it can be used to manufacture 
breaking junctions, where a molecule is pulled out of a 
surface, so that electronic transport measurements can be 
performed^'^. Indeed, the level of geometric control pro- 
vided by the STM even allows mechanical gating experi- 
ments to be carried out, so that a range of different trans- 
port regimes can be investigated for the same molecule^i^. 
In general, both in its topographic and spectroscopical 
modes, an STM measurement consists of a collection of 
different I-V curves for different tip-to-sample positions. 

In the vast majority of cases, calculations of STM cur- 
rents are based on Bardeen's perturbative approach to 
tunneling^ and consist of evaluating the tunneling ma- 
trix elements between the STM tip and the sample. A 
simplified form due to Tersoff and Hamann (TH)^ which 
reduces to the calculation of the local density of states, 
is now mainstream and usually implemented in standard 
electronic structure codes. The TH scheme assumes that 
the STM tip is sufficiently far from the molecule so as 
not to affect its electronic structure. As a result of this 
assumption, there is no self-consistent evaluation of the 
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potential drop nor of the charge re-distribution between 
the tip and the sample. These methods are thus not 
reliable when the tip is close to the molecule and devi- 
ation from the TH current are expected. It is therefore 
important to explore, for STM simulations, the use of 
self-consistent transport methodsi^iiiii^ such as the non- 
equilibrium Green's function formalism (NEGF)^^ com- 
bined with ah initio electronic structures obtained from 
density functional theory (DFT)i^ii^. This is the goal 
of our work, which is based on the use of the electronic 
transport code SMEAGOlJ^iii. 

It is important to realize that NEGF-DFT is com- 
pletely complementary to the standard TH scheme. In- 
deed, it presents problems in the limit of large tip- to- 
sample separation, where TH is most effective. The main 
reason for this is rooted in the fact that SMEAGOL em- 
ploys the numerical implementation of DFT contained in 
the SIESTA code^. The siesta basis set is constructed 
from numerical orbitals whose radial component is trun- 
cated beyond a certain cutoff radiu^, and extended vac- 
uum regions are usually only poorly described. In fact 
in the extreme case of a vacuum region extending well 
beyond all the cutoff radii, all matrix elements between 
the tip and the sample will vanish and the tunneling cur- 
rent will become identically zero. Improvements can be 
obtained by populating the vacuum region with empty 
orbitals (basis functions not associated to a pseudopo- 
tential). However for sensitive calculations of tiny tun- 
neling currents the actual position of the empty orbitals 
can deeply affect the results, most typically by creating 
spurious current oscillations as a function of the tip-to- 
sample distance. Therefore we limit our method to tip- 
to-sample separation smaller than 6A for which reliable 
calculations do not need empty orbitals in the vacuum. 
We call this limit the quasi-contact limit. 

An example of experiments where the tip-molecule in- 
teractions are believed to be important, is that performed 
by Pflaum et. alJ^, in which a monolayer of alkanethiol 
molecules is deposited on a gold surface, and the trans- 
port properties are then probed by using an STM setup. 
The obtained zero-bias conductance is low, being of the 
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FIG. 1: (Color on line) Decanethiol molecule with (a) CH3 
(b) CF3 endgroup. Color code: C=black, S=brown, H=blue, 
F=purple. 

order of 10~^Go at zero bias (the quantum conductance 
Go is defined as 2e^//i, with e the electron charge and 
h the Planck constant). The I-V curves are asymmet- 
ric and this asymmetry increases noticeably when the 
terminating CH3 group is replaced by a CF3 one. It is 
speculated that the additional asymmetry of the CF3- 
terminated molecule is due to a rearrangement of the 
charge distribution near the end of the molecule caused 
by the high electronegativity of the F atoms. This in 
turn generates an electrostatic force between the STM 
tip and the molecule, the direction of which depends on 
the electrostatic potential at the tip. Such interaction 
causes the molecule to be either repelled or attracted by 
the tip depending on the bias polarity, and creates the 
asymmetry in the I-V. 

Calculations using the TH method have been per- 
formed for pentanethiol molecules on gold^^. However, 
the nature of the TH method requires a large tip-molecule 
separation, and the current obtained in these calculations 
is an order of magnitude smaller than that observed in 
the experiments. In this work, the mechanism behind 
the asymmetry in the I-V curves of CH3 and CF3 ter- 
minated alkanethiol molecules on gold is explored with 
NEGF-DFT. The issue related to the electrostatic inter- 
action between the tip and the sample is investigated in 
an approximate way. First we calculate from equilibrium 
DFT the electrical dipole on the molecules. Then we use 
classical electrostatic theory to determine the equilibrium 
position of the molecule with respect to the tip. Finally 
new transport calculations are carried out for the newly 
estimated atomic positions. 



II. METHOD: DFT-BASED 
NON-EQUILIBRIUM GREEN'S FUNCTION 
FORMALISM 

STM transport measurements are equivalent to calcu- 
lating the two-probe I-V curve of a molecule sandwiched 
in between two metallic electrodes. One electrode repre- 
sents the surface to which the molecule is attached and 
the other is the STM tip. The NEGF scheme partitions 
such system into three regions, respectively the two cur- 
rent/voltage electrodes (leads) and a middle region called 



the scattering region (SR) (see figure [2]). 
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FIG. 2: Schematic diagram of a metal- molecule- metal junc- 
tion. A scattering region (the part of the system enclosed by 
the dashed box) is sandwiched between two current /volt age 
probes kept at the chemical potentials /ii and /i2 respectively. 
The electrodes are modeled as being periodic in the direction 
of transport. A number of layers of the electrodes are in- 
cluded in the scattering region to allow the charge density to 
converge to the bulk value. 

The SR includes both the molecule and a portion of the 
leads, and it is extended enough so that charge density 
calculated at the most external atomic layer resembles 
that of the bulk electrodes. The leads, which are assumed 
to be periodic crystals in the transport direction, are kept 
at different chemical potentials /ii/2 = ± eV/2 where 
V is the applied potential bias. We note this is some- 
what different from STM experiments, where typically 
the chemical potential of the surface remains fixed while 
a potential bias is applied to the probe. However, such 
a difference in the definition does not effect the shape of 
the I-V curve, since the final potential drop is calculated 
self-consistently. The SR is described by a Hamiltonian 
Hq. This is used to construct the non-equilibrium Green's 
function 

G{E) = lirn UE + iri) - ff, - Ei - Ss]-^ , (1) 

77^0 

where are the self-energies for the leads, constructed 
by semi-analytical derivation^^ from the electronic struc- 
ture of the bulk. G{E) enters in a self-consistent proce- 
dure to calculate the density matrix, p, of the SR, and 
hence the two-probe current of the deyicei^^^^^. 
The Hamiltonian for the scattering region is generally 
assumed to be a function of the non-equilibrium charge 
density p, which is calculated following the NEGF pre- 
scription as 

1 P"^ 

P=^J G{E)[Trf{E,fir)^T2f{E,fi2)]G\E)dE, 

(2) 

where T1/2 = i[^i/2 — ^1/2]- practice, this integral is 
performed by splitting it into two pdiit^^^^^: an 
equilibrium part to be performed along a contour in the 
complex plane, and a non-equilibrium part to be evalu- 
ated along the real energy axis but that contributes only 
around E^. Finally, the converged Green's function can 
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be used to calculate the two-probe current / through the 
device 
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(3) 

This is effectively the integral between /ii and /i2 (the 
bias window) of the transmission coefficients T{E) = 
Tr[G(^)riGt(E)r2]. 

The NEGF scheme is general and not related to a 
particular electronic structure theory. In the case of 
SMEAGOlJ^iii, the electronic structure method used is 
density functional theory (DFT)^^ in the SIESTA imple- 
mentation. In particular in this work we will use the local 
density approximationi^ (LDA) of the exchange and cor- 
relation functional. 



The orbital resolved density of states (DOS) for iso- 
lated decanethiol both CII3- and CFa-terminated is 
shown in figure [H The energy gap between the high- 
est occupied molecular orbital (HOMO) and the lowest 
unoccupied molecular orbital (LUMO) is in both cases 
rather large (~ 5 eV )^. Most importantly the DOS of 
the two frontier molecular orbitals do not have large am- 
plitude over the C atoms forming the decanethiol, but are 
rather localized on the tt states of S perpendicular to the 
molecule axis (see figures [3fb) and[3l^d)). In the case of 
CFa-termination the F contribution to the DOS is only 
confined to part of the LUMO and to low energy HOMO 
levels. The orbital nature of the HOMO and LUMO is 
further investigated in figure S] where we present the lo- 
cal DOS for the CHa-terminated molecule (the one for 
the CFs-terminated case is rather similar and it is not 
presented here). 



III. RESULTS 

A. Electronic Structure of the Molecules 

The molecules used in the experiments of Pflaum et. 
al. are both methyl-terminated and fluorine-terminated 
alkanethiols. A range of different alkane chain lengths 
were used in that study, but most of the data were col- 
lected for decanethiols (see figure [TJ a) andlljb)). De- 
canethiol consists of an alkane chain of ten carbon atoms 
terminated with thiol (-SH) group. The sulphur atom in 
the thiol group forms a strong bond with gold and an- 
chors the molecule to the surface forming a well ordered 
self- assembled monolayer. 
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FIG. 3: (Color on line) Orbital resolved DOS for the isolated 
decanethiol molecule. Panels (a) and (b) are for the CH3 
termination, while panels (c) and (d) are for the CF3. In 
panel (b) and (d) only the DOS corresponding to tt orbitals 
perpendicular to the molecule axis are displayed. 




FIG. 4: (Color on line) Local density of states for the isolated 
decanethiol molecule terminated with the CH3 group showing: 
(a) the HOMO-1 state (b) the HOMO and (c) the LUMO. The 
local DOS for the molecule terminated with the CF3 group 
is similar and it is not displayed here. Color code: C=black, 
S=brown, H=blue, F=purple. 

We note that both HOMO and LUMO have amplitude 
mainly around the thiol group with little charge spread 
over the alkane chain. One should therefore expect little 
conductance through those states. In contrast the first 
state below the HOMO (HOMO-1) (figure (Ha)) appears 
as a delocalized tt state and it is expected to conduct 
efficiently. 



B. Transport Properties of the Molecules on Gold 
Surface 

In order to calculate the transport properties, we at- 
tach the molecules to the Au fee (111) hollow site via a 
terminating thiol group at a sulphur-surface distance of 
1.9 A^^. This is the equilibrium distance measured for 
the hollow site configuration^^. The arrangement of the 
molecule on the surface and the relative position of the 
tip is shown in figure [5] for the CFa-terminated molecule. 
The configuration for the CHs-terminated case is similar 
and it is not displayed here. In the experiments from 
Pflaum et. al. CHa-terminated decanethiols are tilted 
at an angle of ~ 32° from the vertical direction, while 
for CFa-terminated the tilting angle is ~ 39°—. Our 
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simulation cells however have the molecules placed per- 
pendicular to the surface. This allows the unit cell used 
to be smaller than that required for simulating such tilt- 
ing angles, and the computational overheads are greatly 
reduced. We also make some approximations in mod- 
eling the STM tip. First of all we consider a gold tip 
instead of the used tungsten or platinum-iridium ones, 
in such a way to avoid issues related to having two elec- 
trodes with different work functions. In doing this we 
describe the electronic structure of the electrodes with a 
6s only basis set, which was proved to provide a reason- 
ably accurate choice for transport calculations at only a 
moderate computational cost^*". Secondly we do not con- 
sider an atomically sharp tip but a slightly blunter one. 
This geometry improves the stability of the calculations, 
allowing the self-consistent convergence to be quicker and 
more reliable. 




FIG. 5: (Color on line) Simulation cell used in the transport 
calculations: a CFs-terminated decanethiol molecule is at- 
tached to gold surface and probed with a Au STM tip. The 
case of the Clls-terminated molecule is similar and it is not 
displayed here. Color code: Au=yellow, C=black, S=brown, 
H=blue, F=purple. 

As expected the tunneling current is extremely sensi- 
tive to the distance between the tip and the molecule, so 
that it is rather difficult to determine from experiments 
the precise tip vertical position. Our working strategy 
is then to adjust this distance to match the magnitude 
of the current obtained in experiments. Such a a tuning 
exercise is presented in figure [6] where we show the I- 
V curves for different tip-to- molecule separations, where 
these are defined as the distance between the terminat- 
ing surface of the tip and the C atom in the CH3 group. 
From the figure one can deduce that the best match is 
obtained for a distance of 5. 25 A. The results presented 
in the rest of this section are obtained for this tip-to- 
molecule separation. 

The orbital resolved DOS for both CH3- and CF3- 
terminated decanethiol on Au are shown in figure In 
general and for both the molecule the HOMO is relatively 
close to the gold E^. However, a closer look at its local 
DOS (this is calculated as the local DOS over an energy 
window corresponding to the HOMO level) depicted in 
figures [8] and [9l reveals that such a state is mainly lo- 
calized around the S atoms of the thiol groups. The de- 
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FIG. 6: (Color on line) I-V curves for CHs-terminated de- 
canethiol attached to a gold surface for different distances 
between the C atom in the CH3 group and the probing tip. 
Changing this distance by 0.1 A causes the size of the current 
to change by approximately one order of magnitude at 1 V. 
Panel (b) is a zoom of panel (a) in the current range between 
[-30,20] pA. 




FIG. 7: (Color on line) Orbital resolved DOS for the de- 
canethiol molecule attached to the gold surface. Panels (a) 
and (b) are for the CH3 termination, while panels (c) and (d) 
are for the CF3. In panel (b) and (d) only the DOS corre- 
sponding to TT orbitals perpendicular to the molecule axis are 
displayed for C and S, while in panel (d) all the fluorine tt or- 
bitals are shown. The HOMO-LUMO gap is large, although 
the HOMO is about 1 eV below ^f- 



localized 7r-orbitals capable of efficient charge transport 
across the molecule are further away from the Fermi level 
(approximately 4 eV for both molecules). It is also inter- 
esting to observe the effects that the different terminating 
groups (CH3 and CF3) have over the electronic structure 
of the HOMO-1 level. The most notable difference is as- 
sociated to the charge density distribution along the CH2 
groups forming the decanethiol. While in the case of CH3 
termination the electron charge distributes almost uni- 
formly over CH2 up to near the terminating CH3 group. 
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in the case of CF3 -terminated molecule there is a strong 
modulation with a decay of the charge density as the 
CF3 group is approached. This is due to the Coulomb 
repulsion from the more electronegative CF3 group. Such 
an electrostatic driven charge decay is reminiscent of the 
same effect investigated by STM for molecules on sur- 
faces in the presence of point charges^. 




FIG. 8: (Color on line) Local DOS isosurface for the CH3- 
terminated decanethiol molecule on the Au (111) surface: (a) 
the HOMO-1 level, (b) the HOMO level and (c) the LUMO 
level. Note how the HOMO is localized around the thiol end- 
group and it is not expected to conduct efficiently. Color code: 
Au=yellow, C=black, S=brown, H=blue. 




FIG. 9: (Color on line) Local DOS isosurface for the CF3- 
terminated decanethiol molecule on the Au (111) surface: (a) 
the HOMO-1 level, (b) the HOMO level and (c) the LUMO 
level. Note how the HOMO is localized around the thiol end- 
group and it is not expected to conduct efficiently. Color code: 
Au=yellow, C=black, S=brown, H=blue., F=light blue. 

Despite these differences the absence of any extended 
molecular orbital near the Fermi level of gold is com- 
mon to both molecules and therefore we expect just a 
tiny tunnel current at low bias. This can be appreciated 
by looking at both the transmission coefficients, T{E) 
(figure [TQ|) and the I-V curves (figure [TT]). There is a 
large gap in the resonances in the zero bias transmis- 
sion coefficients of about 5 eV on each side of the Fermi 
level with the edges of such a gap corresponding to the 
HOMO-1 and the LUMO, respectively below and above 
the Au Fermi energy. The transmission of the HOMO 
provides only a small shoulder in T(E) (T ^ 10~^) at 
about 2.5 eV below £^f, confirming that the HOMO does 



not provide an efficient transport channel at resonance. 
However, since the HOMO is the first molecular level to 
enter the bias window (at least for one of the two current 
polarities), it is expected to dominate the low-bias region 
of the I-V. 




FIG. 10: Transmission coefficients at zero bias for decanethiol 
attached to gold surface and terminated with (a) CH3 and (b) 
CF3 endgroup. Note the gap in the transmission of about 5 eV 
on either side of the Fermi level. 

Moving our attention to the I-V curves and the dif- 
ferential conductances shown in figure [Til the most no- 
ticeable feature is their asymmetry. In particular con- 
ductance at positive bias is about between 2 to 3 times 
smaller than that for negative bias. Moreover we ob- 
serve that the conductance of CHa-terminated molecules 
is both larger and more asymmetric of that for CF3- 
terminated. 
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FIG. 11: (Color on line) (a) I-V curve and (b) differential con- 
ductance for decanethiol molecules terminated with both CH3 
and CF3 groups. Note the bias asymmetry, with the conduc- 
tance at negative bias being 2-3 times larger than that for pos- 
itive bias. The conductance is lower for the CFa-terminated 
molecule than that of the CHs-terminated. 

The reason for the conductance asymmetry is rooted 
in the different electronic coupling strength between the 
molecule and the gold surface as compared to the much 
weaker one between the molecule and the STM tip. The 
S atom in fact forms quite a strong bond with the gold 
surface, whereas the tip is separated by a vacuum region 
which provides a substantial potential barrier. Such dif- 
ference in coupling persists to rather close tip-to-molecule 
separations since both the CH3 and CF3 endgroups do 
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not bind to the Au tip. The mechanism leading to 
the asymmetry is then ihustrated in the cartoon of fig- 
ure [T2](a). Since the different coupling strengths of the 
two electrodes with the molecule, electrons are easily 
transferred from the substrate to the molecule, but they 
can hardly hop to the tip. This means that the occupa- 
tion of the molecule is determined by the chemical po- 
tential of the substrate and that the molecule gains elec- 
tronic charges when the electron flux is from the surface 
to the tip (negative bias according to our definition). In 
contrast the molecules lose charge when the flux is in the 
other direction (positive bias). Since charging produces 
a shift of the molecular electronic states towards higher 
energies one expects the HOMO to move into the bias 
window at negative voltages and away from it at positive 
ones. The same behaviour is expected for the HOMO-1, 
while the LUMO follows exactly the opposite pattern (it 
moves into the bias window for positive bias and away 
from it for negative) 



Surface 
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FIG. 12: (a) Schematic energy level diagram of the STM ge- 
ometry. Since the molecule is more strongly coupled to the 
surface than to the tip its occupation is determined by the 
surface chemical potential, i.e. it will charge for negative bias 
and de- charge for positive. Such a charging mechanism pro- 
duces a shift of the occupied molecular level towards higher 
(lower) energy for negative (positive) bias, (b) Cartoon show- 
ing the electrostatic interaction between tip and molecule CF3 
endgroup. In this case there is an electric dipole forming, so 
that the molecule will be attracted or repulsed by the STM 
tip, depending on the bias polarity. 

The electrostatic evolution of the system under bias is 
illustrated in figures [T3l and [Ml where we show the trans- 
mission coefficients as a function of energy for different 
voltages. Clearly the behaviour is the one expected by 
our simple electrostatic model with an upshift (down- 
shift) in energy of the molecular levels at negative (pos- 
itive) bias. Since the highly conducting HOMO-1 and 
LUMO levels never enters the bias window in this volt- 
age range (limited to ±2 Volt), the onset of the current 
is solely determined by the position of the HOMO. This 
however can enter the bias window only for negative bias 



explaining the current asymmetry. Furthermore, since 
the transmission of the HOMO for CHa-terminated de- 
canethiol is larger than that of the CFa-terminated case, 
the conductance of CHa-decanethiol is both larger and 
more asymmetric than that of its CF3 terminated coun- 
terpart. 
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FIG. 13: Transmission coefficients as a function of energy for 
different biases for CHs-terminated decanethiol. Note how 
the resonances in the transmission coefficients due to the oc- 
cupied states move up in energy closer to the bias window 
at negative bias, and move away from the bias window for 
positive bias. The vertical lines mark the bias window. 

Finally, in figure [15] we present the total Miilliken pop- 
ulation, TVm, of the molecule as a function of voltage. The 
figure provides quantitative evidence for the charging of 
the molecule as a function of bias and one can observe 
a total charge variation of about 0.05 e over the entire 
bias range of 4 Volt. Interestingly, and despite the differ- 
ent terminating groups, the charge variation with bias is 
essentially identical for the two terminations. 



C. Influence of the tip-to-molecule electrostatic 
interaction over the transport 

As concluded in the previous section, the asymmetry of 
the conductance with respect to the bias direction found 
in our calculations is due to the asymmetry in the elec- 
tronic coupling of the molecule to the leads. Interestingly 
a strong asymmetry was also observed in the experimen- 
tal conductance measurements for the same moleculesii^. 
However, in the experiments, CFa-terminated molecules 
displayed a far more pronounced asymmetry than the 
CHa-terminated ones. This is in contrast with our find- 
ings where the asymmetry is similar, and in fact it is 
more pronounced for the CHs-terminated case. In agree- 
ment with the initial suggestion from Pflaum et al., here 
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FIG. 14: Transmission coefficients as a function of energy for 
different biases for CFs-terminated decanetfiiol. Note fiow tfie 
resonances in tfie transmission coefficients due to tfie occupied 
states move up in energy closer to tfie bias window at negative 
bias and move away from tfie bias window for positive bias. 
Tfie verticai fines mark tfie bias window. 
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ent magnitude, witfi tfiat associated to CF3 being ttie 
fargest. Finaffy it is important to observe tfiat by farge 
tfie cfiarge distribution over tfie mofecufe is not affected 
by tfie externaf efectric fiefd so tfiat A7V is afmost bias 
independent for botfi tfie terminations. 
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FIG. 16: (Color on fine) Net occupation of the atoms as the 
function of position afong the axis of the mofecufe. 



FIG. 15: Muffiken popufations for decanethiof mofecufe ter- 
minate with (a) CH3 and (b) CF3 groups and attached to the 
gofd surface. Note how in both cases the occupation of the 
mofecufe drops as the bias increases from negative to positive. 



we argue that this additionaf asymmetry is not of efec- 
tronic origin but rather is connected to the efectrostatic 
interaction between the tip and the mofecufe under bias. 
This re-positions the mofecufe with respect to the tip, 
effectivefy changing the tunnefing current. 

The main idea is iffustrated schematicaffy in figure 
[T2lf b). The farge efectronegativity of the F atoms in 
the CF3 group is expected to dispface a substantiaf net 
charge towards the end of the mofecufe. This then inter- 
acts efectrostaticaffy with the surface charge of the STM 
tip, either repeffing or attracting the mofecufe depending 
on the bias pofarity. Such a charge dispfacement is iffus- 
trated in in figure [16] where we show the charge excess 
with respect to neutrafity, A A/", as cafcufated from the 
Muffiken popufation as a function of the position afong 



We estimate the forces acting on the mofecufe by using 
simpfe cfassicaf efectrostatic theory. In the case of CH3- 
terminated decanethiof we cafcufate a force of 3 x 10~^ 
eV in the direction of the tip at the positive bias 
of 2 Voft, and 2 x 10~^ eV in the direction of the 
substrate at the negative bias of 2 Voft. These have to be 
put in refation with the efastic forces needed to stretch 
the mofecufe. DFT cafcufations performed using siesta 
demonstrate that the tip-to-mofecufe efectrostatic force 
at positive bias causes the mofecufe to stretch by about 
0.002 A, whife at negative bias the compression is of the 
order of 0.001 A (at ±2 Voft). Simifarfy, the force acting 
on the CF3 group is cafcufated to be of the order of 3 x 
10~^ eV A~^ in the direction of the substrate when a 
positive bias of 2 Voft is appfied, and of the order of 
4 X 10~^ eV A~^ in the direction of the tip at a negative 
bias of 2 Voft. These forces are in the opposite direction 
to those on the CH3 termination due to the opposite 
directions of the two dipofes. Again DFT cafcufations for 
CF3-decanethiof indicate that the force at positive bias 
causes the mofecufe to compress by about 0.001 A, and 
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the force at negative bias to stretch by 0.002 A. These 
length changes are of the same magnitude for both types 
of endgroup, but they occur at opposite bias directions 
due to the opposite direction of the electrostatic force in 
the two cases. 

The stretching (or compressing) of the molecule un- 
der bias alters the tip-to-molecule separation, which in 
turn affects the magnitude of the current. Let us look 
at the CH3 case first. At positive bias the molecule is 
stretched so that its end is closer to the STM tip, the 
tip-to-molecule separation is reduced and the current will 
increase. In contrast at negative bias, the molecule is 
compressed and a reduction of the current is expected. 
Since the current obtained from our calculations at neg- 
ative bias is larger than that for positive bias, such a 
change in molecule length would contribute to reducing 
the I-V asymmetry in agreement with the experimental 
results. However, a quantitative estimate of the changes 
in the I-V curve as a function of the molecule compres- 
sion/expansion reveals that molecular displacements of 
this tiny magnitude (a few % of an A) are insufficient to 
make a noticeable change to the I-V. 

The results for CFa-terminated decanethiol are simi- 
lar. At positive bias, the molecule is compressed so that 
the CF3 group is further away from the STM tip, while 
at negative bias the molecule is stretched and the tip-to- 
molecule separation is reduced. Since the current calcu- 
lated at negative bias is larger than that at positive bias, 
the electrostatic interaction in this case has the effect 
of enhancing the asymmetry of the I-V curves. This is 
also consistent with the experimental results, showing a 
larger asymmetry in the transport for the molecule termi- 
nated with the CF3 group. However, similar to the CH3- 
terminated case, the estimated molecule length change 
is insufficient to make a noticeable modification to the 
I-V curve. Note that these negative estimates are all ob- 
tained by comparing the geometry at zero bias with that 
corresponding to the largest displacement (at ±2 Volt), 
so that they should be considered as the upper bound 
for the current changes. Therefore we can safely con- 
clude that changes in molecule length can hardly be at 
the origin of the experimentally observed asymmetry in 
the I-V curve. 

A further possibility is to consider changes in the bond- 
ing angle between the molecule and the substrate. As 
mentioned before, decanethiols do not arrange vertically 
to the surface, but they rather form angles of 32° and 
39°, respectively for the CH3 and CF3 termination^. 
Hence it is likely that the interaction between the STM 
tip and the endgroups causes the molecule to rotate ei- 
ther away from or towards the STM tip, as shown in 
figure [171 The direction of the rotation depends again on 
the sign of the electrostatic force. However since the en- 
ergy required to change the bonding angle is expected to 
be rather smaller than that required for an axial distor- 
tion, we expect changes in bond angle to produce larger 
changes in the tip-to-molecule separation. 

Using the same electrostatic forces calculated before we 
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FIG. 17: Schematic illustration showing the tilting angle of 
the molecule over the Au surface and the rotation due to the 
interaction with the STM tip. 



estimate for the CH3 termination a 1° rotation towards 
the STM tip at positive bias (2 Volt), and a 0.5° rotation 
towards the surface at negative bias (-2 Volt). Similarly 
for the CF3 termination, the rotations are along the op- 
posite direction with changes in the bond angle of 0.5° 
towards the surface at positive bias and 1° towards the 
STM tip at negative bias. These rotations are sufficient 
to alter the I-V characteristics, as shown for both types 
of molecule in figure [TSl 
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FIG. 18: (Color on line) I-V curves and differential conduc- 
tance as a function of bias obtained both with and without 
considering the effects of the electrostatic-induced rotation of 
the molecule. Panels (a) and (b) show the I-V and d//dy 
curves for the CH3 termination, and (c) and (d) show the 
same for the CF3 case. 

In the case of CH3 termination the calculated rotation 
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is actually sufficient to suppress the asymmetry in the 
I-V present because of the asymmetric coupling to the 
electrodes. The current at 1 Volt is now larger for posi- 
tive bias than that at the same negative bias. In contrast 
for CFa-terminated molecules, the rotation is such that 
the I-V asymmetry already present is enhanced. Thus, 
when the electrostatic-induced rotation of the molecule is 
taken into account, the I-V curve for the CFa-terminated 
molecule appears as more pronounced than that obtained 
for the CH3 -terminated one, in semi-quantitative agree- 
ment with experiments^^. 

In concluding we wish to point out a few possible 
sources of error in our calculations, which might affect a 
fully quantitative comparison with experiments. Firstly 
we remark that the Au atoms of both the tip and the sub- 
strate are described at the 6s-only level. This provides 
a good description of the Au Fermi surface, thus of the 
transport properties, but usually it does not provide ac- 
curate interatomic forces. Secondly the calculated forces 
between the molecule and the tip certainly depend on 
the actual shape and detailed position of the tip. These 
are difficult to characterize with certainty and errors are 
therefore difficult to avoid. Finally the actual value of the 
tilting angle depends on the inter-molecular interaction 
is the self-assembled mono-layer, and thus the force re- 
quired to change the angle would depend on the details 
of the molecular coverage (density, symmetry, bonding 
site, etc.). 

IV. CONCLUSION 

We have demonstrated that the DFT-NEGF code 
SMEAGOL can be used to simulate STM-type experiments 
in the near to contact regime. This is under the working 
condition of the tip-to-sample distance being sufficiently 
small that the basis orbitals have not been artificially cut 
off (i.e. the vacuum region between the tip and the sur- 
face is still well described). SMEAGOL then allows us to 
investigate systems in which the tip interacts with the 



molecule, and to study effects of finite bias on the elec- 
tronic structure of the molecule. 

Calculations for alkanethiol molecules with STM-type 
arrangements show strong asymmetry in the I-V curves, 
which can be explained by the asymmetry in the coupling 
to the two different electrodes (substrate and tip). How- 
ever, the calculated asymmetry is similar for both CH3- 
and CFa-terminated decanethiols, in contrast to the ex- 
perimental measurements, showing a far stronger asym- 
metry for the case of CF3 termination. We have then 
investigated the original suggestion, which attributed 
the asymmetry to small configurational changes of the 
molecule under bias due to electrostatic interaction be- 
tween the tip and the substrate. Our calculations demon- 
strate that changes to molecule length are too small and 
unlikely to have any major impact on the I-V curves, 
however a rotation of the molecule may cause significant 
changes to the current. 

As a final remark we point out that the results pre- 
sented in this work have obtained by using the LDA. 
However, LDA has several shortfalls which can strongly 
affect electron transport calculations. In this case, ap- 
proximate self-interaction corrections such as ASIG2ii2^ 
are unlikely to offer a substantial improvement, since the 
conductance is due to off-resonance tunneling, and would 
not be particularly sensitive to the exact position of the 
molecular orbitals. However, for these molecules, an ac- 
curate calculation of the electric polarisability beyond the 
LDA might be important^^. 
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